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We find a simple expression for critical curves of a binary gravitational lens. On the 
basis of this, we present a parametric representation of such curves. The caustics can also 
be expressed with the same parameterization. The present result is helpful for efficiently 
constructing many templates of light curves due to binary systems, particularly extrasolar 
planets, which cause spikes in the light curves when a source crosses the caustics. 

§1. Introduction 

Since the pioneering work of Schneider and Weii3 (1986), gravitational lensing 
by binary objects has attracted a great deal of interest. It offers an important 
tool also in the search for extrasolar planets, since it is independent of the 

Doppler method, which cannot determine the planet mass because of the unknown 
inclination angle. ''^ In general, critical curves, on which the lensing amplification 
diverges, play an important role in lensing studies: The number of images changes 
at these curves, and giant arcs are also observed (For a review, see Ref. 2)). To this 
time, the curves for binary gravitational lenses have been drawn by solving a lengthy 
equation. This is a fourth-order equation,*)'^) which can be solved analytically only 
with a quite lengthy form of roots for instance by using Mathematica. Thus it has 
been believed that no simple parametric representation can be obtained. 

It seems possible, however, to introduce a parameter using sine and cosine func- 
tions such that we can determine a simple representation of the curves, because 
neither the sine nor cosine function is algebraic. In the complex formalism, by con- 
trast, a parameter is introduced in the beginning, and the fourth-order equation 
is obtained at the end. Here, we should recall that algebra allows only a finite 
number of operations of addition, subtraction, multiplication, division and root ex- 
tractions. However, introducing a parameter through sine and cosine functions 
changes the situation completely, because these functions themselves contain an in- 
finite number of operations. Hence, our result is not algebraically equivalent to the 
previous one. The main purpose of this paper is to obtain a compact expression of 
the critical curves: We re-examine the equation for them and show that it can be 
rewritten in a tractable form. This is a key discovery, allowing us to parameterize 
the critical curves and caustics in a simple way, as shown in the following. 
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§2. Critical curves and caustics for a binary gravitational lens 

2.1. A lens equation 

We consider a binary system located a distance Dj^ from the observer. It con- 
sists of two point masses, Mi and M2, whose positions are denoted by and L2, 
respectively. The displacement vector from object 1 to 2 is written as L = L2 — L-^. 
For later convenience, let us define the Einstein ring radius angle as 



_ 4GMPLS 

where G is the gravitational constant, M is the total mass. Mi + M2, and Dg and 
-^LS denote the distances between the observer and the source and between the lens 
and the source, respectively. In units of the Einstein ring radius angle 6-^, the lens 
equation reads 

where /3 and 9 denote the vectors representing the positions of the source and image, 
respectively, and we defined the mass ratios fi and 1^2 and the angular vectors for 
the positions and separation £1 , £2 and £ as 



M, 



(2-3) 



' Ml + M2 ' 

for z = 1, 2. Clearly, we have the identity 1/1 + 1/2 = 1. Equation (2-2) is a set of two 
real coupled fifth-order equations for (O^, 9y), which is equivalent to a single complex 
quintic equation for 6x + iOy.^"^^'^^^ It has been shown recently that this equation 
can be reduced to a single real quintic equation. 

2.2. Critical curves 

Critical curves are curves in the lens plane, where the Jacobian of the lens 
mapping \d(3/d0\ vanishes. Now, we define the angular separation between the 
image and the lens i {i = 1,2) as 

Pr = \e-£il (2-6) 

where pi does not vanish, because the location of the image is different from those 
of the lens objects. 

Let us denote the height of the image measured from the separation axis by h. 
Then, straightforward computations of a vanishing Jacobian give us 

p\pi = {^2pI + viplf - Mz^2^'^'. (2-7) 
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By rearranging the expression for the separation as ^ = \J pi — hP' + \^ p\ — }iP^ we 
find 

,2 4pf p1 - {e -p\- pif 

^ 4^2 • (2-8) 

Substituting this into Eq. (2-7), we obtain 

p\pi = MpI - ^i^^f + yi{pl - ^2^^?. (2-9) 

where we defined ^= Here, it should be noted that this expression is independent 
of the choice of the origin of our coordinates. In the complex formalism, for instance, 
the equation for the critical curves is inevitably a fourth-order polynomial equation, 
whose coefficients depend on the choice of the origin of the coordinates, so that 
we need a numerical method to solve the quartic, and the resultant expressions are 
quite lengthy. 

Our procedure for obtaining a parametric expression of critical curves is divided 
into two steps. First, we find a simple parametric expression of pi and p2- Next, 
using this expression, we obtain those of 6x and dy. For convenience, we define 

u = -^, (2-10) 



/V2 

SO that Eq. (2-9) is rewritten simply as 



^2 (2-11) 



2 +^ 2-^' (2-12) 

V u 



which allows us to introduce a parameter s G [0, 27r) according to the relations 



1 



coss = (2-13) 

V 



1 _ v^^^ 

sins= (2-14) 



Equation (2-13) can be rewritten as 

1 cos s 



(2-15) 



Substitution of this into Eq. (2-14) gives us a quadratic equation for u: 

(sin s)u'^ - {y/U^f sin s - ^/l^f cos s + l)u + ^/IT^f = 0. (2-16) 
Equation (2-16) is solved as 

a±Vb 

2 sm s 
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where we defined 



cos a = (2-18) 

sin a = (2-19) 

A = l+fsm{s-a), (2-20) 

B = -Afcosasins. (2-21) 



Substituting Eq. (2-17) into u of Eq. (2-13) leads to 

c±Vb 



(2-22) 



2 cos s 
where we defined 

C = l-fsm{s-a). (2-23) 

Up to this point, we have not specified the coordinate system. Now, let us 
introduce polar coordinates (r, 0) whose origin is located at the mass Mi and whose 
angle is measured from the separation axis of the binary. This implies that = 
and £2 = I- In terms of the polar coordinates, p\ and p2 are expressed as 

p\ = r\ (2-24) 
pi = r^-2^rcos(^ + /^ (2-25) 

Then, Eqs. (2-24) and (2-25) are solved for r and cos0 as 



- ,1/4 



u, (2-26) 



'""^^ o^V^A- ' ^^^^ 

where wc used Eqs. (2-10), (2-11), (2-17) and (2-22). Hence, in terms of u and v, 
the critical curves are expressed as 



9x = r cos (j) 

_ \fv\u — ^/viv + 



(2-28) 



±r\/l — cos 



2 , 



where we defined 

t + l^xl^u + 2^fv + (2-30) 

We are now in a position to consider the constraint on pi and that the two 
lens objects and the image form a triangle. This condition is expressed as 



\pi-p2\<t, (2-31) 
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which can be rewritten as 



1/4 



(2-32) 



Since u and v are positive, we can show that this inequahty is equivalent to the 

relation D > 0. 

It is possible in principle to determine allowed regions for s by analytically 
treating Eq. (2-32), which is the condition for a triangle. In practice, however, it is 
much easier to draw the critical curves for s G [0, 27r) only if we ignore s for D < 0. 

In the case of a symmetric binary, i.e. ui = = 1/2, the condition D > 
becomes 

(sin s + cos s)^ + 2 V2(sin s + cos s) - 4 < 0. (2-33) 



This is solved as 



sin s + cos s < 



(2-34) 



where we used the fact that both | sins| and | coss| are less than unity. Hence, we 
find 

^-sin-i(V3-l) <u< ^ + sin-i(V3-l), (2-35) 

where we assumed s G [7r/2, 57r/2) for convenience. Otherwise, for [0, 2tt), the allowed 
regions arc apparently divided into two, s G [0, sin~^(\/3 — 1) — 7r/4) and [37r/4 — 
sin-i(V3- l),27r). 
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Fig. 1. Critical curves for a binary gravitational lens with the separation of the Einstein ring ra- 
dius £ = 1. The dotted and dashed curves correspond to the mass ratios fi = 1/2 and 1/5, 
respectively. 
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Fig. 2. Critical curves for a binary gravitational lens with the separation of the Einstein ring di- 
ameter £ = 2. The dotted and solid curves correspond to the mass ratios ui = 1/2 and 1/5, 
respectively. 



As an illustration, critical curves are given for two cases, that in which the binary 
separation equals the Einstein ring radius (Fig.l), and that in which the separation 
equals the Einstein ring diameter (Fig. 2). These figures demonstrate that the value of 
the binary separation is of crucial importance in determining topological properties 
of the critical curves, in particular, whether there is one loop or two, and that a 
non-half mass ratio causes the curves to be asymmetric. ^^'^^ 

2.3. Caustics 

The parametric representation of critical curves is given by Eqs. (2-28) and 
(2-29). Substitution of these into the lens equation (2-2) with .^i = and t2 = i 
gives immediately a representation of caustics in terms of the parameter s, 

where 0{s) is given by Eqs. (2-28) and (2-29). 

The caustics corresponding to Figs. 1 and 2 are given in Figs. 3 and 4, which 
show that the binary separation is important in determining the topological features 
of caustics, while caustics are distorted due to the mass ratio. 
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Fig. 3. Caustics for a binary gravitational lens with the separation of the Einstein ring radius £ = 1. 
Each curve corresponds to that in Fig. 1. 
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Fig. 4. Caustics for a binary gravitational lens with the separation of the Einstein ring diameter 
i = 2. Each curve corresponds to that in Fig. 2. 
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§3. Conclusion 



First, an equation for the critical curves of a binary gravitational lens was recast 
into the compact form (2-9). On the basis of this, we obtained a parametric represen- 
tation of these curves as Eqs.(2-28) and (2-29). A parametric representation of the 
caustics was also obtained as Eq. (2-36). These representations are helpful to study 
binary Icnsing, particularly for a number of sets of the mass ratio and separation. 

For instance, when we discover extrasolar planets by microlensing, a standard 
strategy is to detect a slight change in light curves using a small telescope, and 
then to follow this up using a large telescope on the ground or in space as soon 
as possible. In order to be able to respond as quickly as possible, it is necessary 
to prepare many light curve templates, in which spikes appear owing to planetary 
microlensing. The positions of such spikes, which correspond to crossing caustics, are 
essential for templates used in initial observations. The detailed shapes of the spikes 
are important to allow for analysis of light curves with follow-up observations. The 
expression obtained in this paper may be helpful for constructing these templates 
with good accuracy and efficiency. This is a future subject. 
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